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Abstract 

We discuss the correlations of charming penguin contributions to ^ vrvr, t:K and KK using 
approximate flavor SU (3) symmetry. Strong constraints are found from the direct CP asymmetries 
especially in ttK modes. We make a global fit to the latest data and find that only a small charming 
penguin is allowed, and the size of color-suppressed tree amplitude (C) relative to tree amplitudes 
(T) remains large C/T ~ 0.6, which disfavors the possibility of a large charming penguin as an 
explanation for the vrvr puzzle. We show that a small charming penguin can still have sizable effect 
in the time-dependence CP asymmetries in KK mode. 
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I. INTRODUCTION 



With the successful running of the two B factories, the B physics has entered a precision 
era. Although the current data of hadronic B decays has show an overall agreement with 
the Standard Model (SM), there are a number of modes with unexpected decay rates and 
CP asymmetries, which are often referred to as puzzles. One of the puzzles in vrvr modes is 
a large averaged branching ratio of vr^vr^ relative to vr+vr", the current data read [l| 

= D ^ = 0-51 ± 0.08 1 

which is significantly larger than theoretical estimations. Another one in nK is the difference 
in two direct CP asymmetries 

A^p{7^+K-) = -0.097 ± 0.012 , Acp{n^K~) = 0.050 ± 0.025, (2) 

Both of the puzzles require a large color-suppressed tree amplitudes in flavor SU(3) topology, 
which is difficult to obtain from short-distance contributions. So far a satisfactory explana- 
tion is not yet available. There are other potential puzzles regarding the branching ratios 
and time-dependent CP asymmetries in nK mode which are relevant to the possibility of 
new physics. In the present work, we focus on the former ones which are more relevant to 
the hadronic dynamics. 

ong-distance Final State Interactions (FSIs) may 



It was emphasized in literature that 
play important roles in these modes [2, S, 4\, such as the inelastic rescattering chan- 
nel B -D-D(s) -mi{K), KK at meson level. Topologically they are equivalent to 
the charm-quark loops in the contractions of local operators Ql = {dc)v~A{cb)v-A8i'nd 



m. 



Q2 — {daCf5)v-A{cf^ba)v-A at quark level which are referred to as charming penguins[ 
Experimentally Br{B D+D-^^) = [1.9 ± 0.6(65 ± 21)] x lO-^fl are about 40(300) times 
larger than that of B ^ n^n^^K^). As a consequence, a tiny OZI violating DD — > vrvr 
ead to significant changes in branching ratios and CP asymmetries in iriT^nK) 
llll |. A large charming penguin with an appropriate strong phase may simul- 
taneously suppress Br{n~^'7i~) while enhance _Br(7r°7r°), thus providing a solution to the tttt 
puzzle. 

The effects of charming loop have been discussed at both quark level and meson level. The 



mixing inav 



modes p, llO 



situation is not yet conclusive. Estimations based on pQCD 



12| and QCD sum rules [l3| favor 
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a small size. While in the framework of Soft Collinear Effective Theory (SCET) the charming 
penguin could be large, depending on the jet functionjl^. The meson level calculations using 
effective Lagrangian for mesons favor a large charming penguin comparable to QCD penguin 



in ttK 
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16|, 117|, Q . But the patterns in vrvr data can not be well reproduced 
Note that a global analysis using approximate flavor SU(3) symmetry for all the nn, nK 
and KK modes may provide a powerful constraint on charming penguins. This is because 
the presence of charming penguin not only modifies the individual decay amplitude but also 
changes the correlations among them, which has not been enough emphasized in previous 
analysis. The correlations are of particular importance in distinguishing charming penguin 
from QCD penguin. First, although the two type of amplitudes always appear together, in 
AS = 1 modes they are both nearly real, but in AS = modes they differ by a phase angle 
/3 of the unitarity triangle. The correlations in the predictions of direct CP asymmetries are 
changed. Second, the AS = 1, B ^ nK modes are penguin-dominant, which constrain the 
absolute size of charming penguin together with QCD penguin, while the AS = 0, B ^ mi 
modes are tree-dominant and more sensitive to the tree-penguin interference. A strong 
constraint comes when they are combined together. Finally, the AS = 0, 5 — KK modes 
provides a testing ground for the charming penguin. In the SM, the time-dependent CP 
asymmetry S{KsKs) is nearly zero because only QCD penguin contributes. The presence of 
charming penguin provides an additional amplitude with different weak and strong phases. 
Thus a significant deviation from zero is possible. 

The present work is organized as follows. In section [Tll we discuss the nontrivial corre- 
lations caused by long-distance charming penguin using the QCD factorization results for 
short-distance contributions. In section IIIIl we make a largely model-independent global 
determination for the charming penguin using the latest data. The results show that a 
small charming penguin is favored, which can not play any significant role in resolving the 
TTTT puzzle. But it can still significantly affect the prediction for S{KsKs)- Some remarks 
and conclusions are in section 



II. CHARMING PENGUIN CONTRIBUTION TO INDIVIDUAL MODES 

The simplest way to see the correlated contributions from charming penguins in different 
modes is to fix other hadronic amplitudes to their theoretical values. To this end, we 
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decompose the whole decay amplitudes (A) into short-distance (Asd) and long-distance 
{•^sl) part A = AsD + Ald, and take the short-distance part from theoretical calculations. 
The long-distance part is assumed to be dominated by charming penguins. The decay 
amplitudes are related to the observable of decay branching ratio and direct CP asymmetry 
as follows 

where Pc is the momentum of final state meson in the B meson rest frame and tb = 
1.530(1.638) X lO^^^s [sl is the neutral (charged) B meson life-time. The time dependent 
CP asymmetry is 

r(50 ^ fcp) - T{B^ ^ fcp) 



acp{t) 



T{B<^^fcp) + TiB<^^fcp) 

S ■ sm{AmB -t) - C ■ cos{AmB -t). (4) 



The definition of quantities S and C are given by 

where (q/p) = e"^*'^ in the SM with P one of the anlges of the unitarity triangle (UTV. In 



what follows we take the CKM matrix elements Vub and Vet from the global CKM fits 

Kb = (3.57 ± 0.17) X 10-^ = 0.0405;°:°°^. (6) 
To fix the profile of the UT we also use the best fitted value of {l^ 

7 = l.nOt'ofrl (7) 

which corresponds to a best fitted P = 0.379 ± 0.017. 

Recently the theoretical calculations for hadronic matrix elements have been improved 



'actorization for spectator 
~l 

2j|. In QCD factorization 



to next to leadin g a 't order (NLO) in the framework of QCD 
scaterings and in perturbative QCD (pQCD) 

approach, the hard spectator scattering effects can lift a cancellation between leading term 
and vertex corrections, resulting in a significant enhancement in the effective coefficient 
a2(7rvr) by a factor of ~ 3 and improve the agreement with the data. Nevertheless generating 
a large enough spectator scattering effects still require tuning of input parameters and the 



latest calculation still favor a Br^n^n^) lower than the current data [2l||. It remains to be 



seen if there is futher enhancement from NNLO calculations [25|, l26|. Note that in the pQCD 
approach, although the NLO results improve the predictions for the direct CP asymmetries 
in irK modes, there is no significant enhancement found in vr^vr". 



The whole charmless B deca y a rnp 



quark flavor flow diagrams 
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28|, 
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itudes can be described by a set of flavor topological 



3l| . In this approach the decay amplitudes are 



expressed in terms of diagrams such as tree (T), color-suppressed tree (C), QCD penguin 
(P and Pt„), electroweak penguin (PevkO, color-suppressed electroweak penguin (Pew) ^t*^- 
In the presence of charming penguin ( denoted by D for DD intermediate states ), the decay 
amplitudes for vrvr modes are given by 



-^(vr+TT-) = A„(T + E-Ptu-PA- -Pi 



EW 



\c{P - D 



Pa + -^Pew)-> 



1 



^(tt^tt") = -=[A„(C -E + P^^ + Pa-P, 



1 

7! 



EW 



i;Pew) + UP -D + Pa 



P 



EW 



-P? 



EW)\^ 



)], 



[Au(T + C — Pew — Pew) ~ ^c{Pew + Pew)] 



(8) 



The CKM factors are defined as = V*^(^g-^Vqh. In general, the QCD penguin has three 
part A„Pu + AcPc + A^Pf, which is recombined as Ptu = Pt — Pu and P = Ptc = Pt — Pc- The 
amplitudes T, C and P etc can be calculated and the typical values (in units of iC^eV) from 
QCD factorization are 



T — A^^j^ai 



~ 0.89 - 0.02i, 



0.24- 0.02i, 



Ptu 

p 



0.076 + 0.0292, 
0.084 + 0.015Z, 



(9) 



where A, 



GfIttPq 



[m 



B 



m 



l)/\/2. The numerical values in the above expressions 



are in accordance with the central values of NLO effective coefficients in QCD factorization 
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approach 





= 0.975^°:^?^ + 


(- 


-0.017+°:^-)^, 






= 0.275l°:?i + 


(- 


-0.024^°:^^?)^, 




"4,7r7r 


= -0.024^°:°°^ 


+ 


(-0.012^°::^; 


i 




= -0.0281°:°°^ 


+ 


(-0.006^°:°°^; 


i 




= -0.0601°:°?^ 


+ 


(-0.020^°::^; 


i 


X 6,7r7r 


= -0.065^°:°- 


+ 


(-0.0101°:°°^; 


i 



The short-distance calculations suggest a t— quark dominance in QCD penguin such that 
Ptu — P, and tiny annihilation type diagrams E, A and Pa which are power suppressed. 



A. TTTT modes 

We begin with a re-examination of tttt puzzle in the presence of charming penguin D. In 
the limit of T,C ^ P, D, the ratio Rt,tt can be expanded as follows 

E ~ — 



P P 

1 + 2 (1 - cos 7) ( — cos((5t - Sp) + cos(5c' - 



+ 2ujrD { j; cos((5t - 5d) + ^ cos{Sc - 5d) ) cos 7 



(11) 



where u = \Xc/^u\ — 2.73, and td = D/P is the size of charming penguin relative to QCD 
penguin. It is evident that the charming penguin has opposite contributions to tt+tt^ and 
7r°7r° modes. In order to enhance i^^r^ one needs cos{6c — 5d) > and a large r^j. In Fig. 
1 we plot the ratio Rt,t, as a function of with different strong phases. In the numerical 
calculations we use the full expressions for decay rates and CP asymmetries. 

It is shown in the figure that for a typically small strong phase 6d = 30°, a large 2 < < 
2.5 is needed to meet data of R-„t,. For large strong phase 6d > 90° , an even larger vd > 3 
is required. This confirms previous phenomenological studies in favor of large charming 
penguin. The direct CP asymmetric measurements provide different constraints. In the 
limit T, C 3> P, -D, the direct CP asymmetries are approximated by 

p 

acp(7r+7r ) ~ 2uj— (sin(5T - 5p) - td sin(5T - So)) sin 7, 
p 

acp(7r°7r°) ~ -2c^^ (sin(5c' - ^p) - ro sm{6c - Sd)) sin 7. (12) 
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Figure 1: The ratio R-,^t^ as function of charming penguin. Four curves corresponds to the strong 
phase 5d = 30°(sohd), 60° (dashed), 90° (dotted) and 120° (dot-dashed) respectively. Other param- 
eters are default in QCD factorization estimations. 



Similar to the vrvr decay rates, the charming penguin contributions to the direct CP asym- 
metries are again opposite. Since 5t and 5c are small, roughly speaking for < 5/) < 180°, 
it enhances acp(7r+7r^) while suppresses acp(7r°7r°) to negative values. The numerical results 
are shown in Fig. 2. 

For small phase 5d = 30°, the current data of acp(7r+7r^) restricts the size of ro to be 
1.5 < r£, < 2. Note that there is still a significant difference between two B factories on the 



measurement of acp(7r+7r" 



3d, 



33| 



0.21 ± 0.09 ± 0.02(5a6ar) 
0.55 ± 0.08 ±0.05(Se//e) 



(13) 



The Babar measurement favors a smaller acp(vr"'"7r') and the constraints on rr, is stronger. 
Note that the constraints on the size of td and the strong phase dr, from the acpij^^n^) and 
-R^TT are opposite. The preliminary data although with large uncertainty are in favor of a 
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Figure 2: acp('/r+7r~) and Ocp (vr'^Tr'') as function of D. Four curves corresponds to the strong phase 
5d = 30°(sohd), 60° (dashed), 90° (dotted) and 120° (dot-dashed) respectively. Other parameters 
are default in QCD factorization estimations. 

positive acp(vr''7r°), which disfavor any large value of td with 5d in the range (0,7r). 

A more significant dependence can be seen in the time-dependent CP asymmetries 
which are approximated by 

p 

S{'K^'K ) — sin2(/3 + 7) + 2^7— (cos((5p — 5t) — ro cos{Sd — St)) sin7cos2(/3 + 7), 
5(7rV) ~ -sin2(/5 + 7) -2cu^(cos(5p-5c) -rDCOs(5D-5c))sin7cos2(/5 + 7). (14) 

Since the current global CKM fitting give a /3 + 7 close to 7r/2, the leading term is suppressed 
for both TT^TT^ and which significantly enhances the charming penguin effects. As 

shown in FigjSl for 6d < 60° the data of S'(7r+7r~) exclude the possibility of charm penguin 
since the short distance contribution is already above the experiments. The data of ^(Tr+Tr^) 
favors a larger strong phase 6d. The charming penguin contribution can be either positive 
and negative, depending on cos{6r) — 6t)- Note that in vr^vr^ mode the charming penguin 
contribution is proportional to P/C much larger than that in tt+tt^ which is proportional to 
P/T. Thus ^(Tr^Tr^) has the strongest charming penguin dependence, which can be clearly 
seen from FigjSl For 6d = 30° and ro = 1, the value of S^n^n^) can be reduced to around 
zero. The future precision measurement of S'(7r''7r'^) will provide a very strong constraint on 
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Figure 3: ^(vr+vr^) and S'(7r'^7r'^) as function of rjj. Four curves corresponds to the strong phase 
6d = 30°(sohd), 60° (dashed), 90° (dotted) and 120° (dot-dashed) respectively. Other parameters 
are default in QCD factorization estimations. 

B. ttK modes 

We process to discuss the nK modes. The diagrammatic amplitudes for ttK modes are 
given by 

-A{n+K~) = KiT - Ptu - \pew) - K{P -D + \p^w) 
-^(7r°i?°) = -^\K{C + - Pew - \p^w) + K{P - D - P^w - \p^w)] 

A{ii-K') = K{A - P,, + Ipf^) - K{P -D- Ipf^) 
-^(7r°i^-) = ^[A«(T + C + A-P,^- Pew - \p^w) ' UP - D + Pew + ^PEw)] 

(15) 

where the amplitudes T, C and P for ttK modes can be obtained by replacing A^,.,, into A^^k- 
The electroweak penguin amplitude is calculated from the effective coefficients aj^T^x and 

Pew = ^A^Kial^K ' «9,.i^) - (1-317 + O.OlSz) x 10-^ 
PEW,tu = ^A,i^(4,,^ - al^j,) ~ (1.315 - 0.014^) x lO'^ (16) 
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which corresponds to the effective coefficients 



2l| 



a^,,;, = 0.058t°;°f7 + (0.015 



-O.OIOn 
-0.006 J 



-0 819+° °^° 



(0.029 



+0.053N 



c 



0.0101^:^1^ + (o.ooo: 



-0.003>, • 
-0.006 J ^ 



-0.023/^ ' '^9,ttK 



+ m^^trors)^- (17) 



One sees again that the t— quark dominance leads to Pew — Pewm- what follows we 
neglect teh subleading color-suppressed diagram Pew- modes one can define the 

following two ratios for neutral and charged modes [jj 

Br(n+K-) „„„,„„_ „ 2Br{TX^K-^ 



Rr, 



0.98 ±0.07, Rr 



1.12 ±0.07. 



(18) 



The penguin dominance leads to an estimation of Rn ^ Rr. ^ 1. Due to a cancelation in the 



subleading term, Rn — Rc holds to a high accuracy [3^, 135||. Although earlier data showed a 
small Rn, which is usually referred to as nK puzzle, the latest measurements show that this 
puzzle has been significantly reduced. It is easy to see that the dominant charming penguin 
contributions cancel out in both Rn and Rc and the remaining parts are suppressed by a 



CKM factor ^ = 
Rn ^ R, 



\K/K 



0.02 

T 



^1-2^ 
C 

-pfD cos{5c - Sd) 



C T 
p cos((5t - (^p) ± -p cos(5c - h) - -pTD cos(St - Sd) 



COS7 ± 2— p- cos(dp^^ - bp) 



(19) 



The previous global fits without charming pen guin show a remarkable agreement between 



theory and experiment in penguin amplitudes 
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3- 



The pure penguin mode such 



as -K'K' constrain strongly the absolute values of \Pe — De^^\. However, the charming 
penguin can still be sizable when it carries large relative strong phase. In this case, due to tiny 
but nonzero weak phase difference between D and P, unacceptably large CP asymmetries 
can be induced when the charming penguin and QCD penguin are comparable in size. The 
charming penguin contributions to the direct CP asymmetries in the limit D -C P are given 

by 



acp{TT^K ) 

acp(vr-i?°) 
acpin'^K') 



T 

-2^— [sin(5T - Sp) - td sin(5r - 5d)] sin 7, 
C 

2^— [sin(5(7 - bp) - td sin(5c - 5d)] sin 7, 
2^rD sin(5z) - 5p) sin 7, 



T 



C 

sm{6T - 6p) + — sm{6c - Sp) - td sm{ST - Sd) 



sm7. 



(20) 



10 



Numerical calculations for the direct CP asymmetries are given in figlH which shows that 
when the value of P is fixed by the QCD factorization default value there is little room 
for D except for the unreasonable region D ^ P. The strongest constraint comes from 
acp(7r+_ft'~), for 6d = 30° the allowed value of ro is very narrow around re ~ 0.7. Large ro 
is only allowed for some special settings such as 6d = 120°. The data of acp(7r°fC^) imposes 
a similar constraint, for 6d = 30° the allowed vd is around 0.5. For other values of the 
strong phase 6d = 30° ~ 120° the allowed td is even smaller around 0.3. For the other two 
modes 7r°X° and n~K^ the constraints are much weaker due to the weakened or vanishing 
tree-penguin interferences. The data only disfavor the value of ~ 1. The correlations 
among the four modes can be clearly seen from figlH The charming penguin contribution 
in acp{7!'^K~) is opposite to those in acp(vr°_ft'°) and acp(vr~K°) but similar to a{7i^K~). For 
small 6d and td ~ 1, acp{7!'~^K~) and a^n^K') reach their minimum, while acp(7r°^°) and 
0'cp{'^~t^^) close to their maximum. 

Note that the charming penguin in AS" = 1 modes are nearly CP conserving, the time- 
dependent CP asymmetry for n^Ks remain unchanged for small td 



This is due to the fact that the charming penguin contribution cancels in the ratio A/ A at 
the leading order. 

C. KK modes 

The decay amplitudes for KK modes are given by 



The K^K^ and K~K^ are pure AS* = penguin modes. The charming penguin contribution 
is similar to that of n'K^ except for the CKM factors. From the latest data, the QCD 
penguin amplitudes can be extracted and found consistent with that from n~K^. Note that 




(21) 



A{K+K') = -\^iE + PA) 

A{K'K') = -K{Ptu - \p^w) - UP -D- \p^y,) 
A{K-K') = UA - Ptu + ipf^) - A,(P -D- ipf^) 



(22) 
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Figure 4: acp{'7T~^K~), acpiir^K^), acpiir^K^) and acp{iT^K~) as functions of r£). Four curves 
corresponds to the strong phase 6d = 30°(sohd), 60° (dashed), 90° (dotted) and 120° (dot-dashed) 
respectively. Other parameters are default in QCD factorization estimations. 

in the SU(3) limit, the direct CP asymmetry of K'^K^ is directly linked to the '7i~K^ 



uS^ — 2uj cos 7 + 1 



acp{K'>K'>) ~ -2— — ^ _ sm{SD - Sp) sin 7 



acp{T^-K'') (23) 



,^(u;^ — 2^7 cos 7 + 1) 

Note that the asymmetry is enhanced by a factor 1/^ ~ 50. From the current \a bound 
acp(vr-i?0) = 0.009 ± 0.025, acp[K^K^) can easily reach to C(-0.5) for ~ 0.5. A 
stronger td dependence can be found in the time-dependent CP asymmetry of KsKs- For 
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Figure 5: acp{K^ K^) and S{KsKs) as function of r/). Four curves corresponds to the strong phase 
6d = 30°(sohd), 60° (dashed), 90° (dotted) and 120° (dot-dashed) respectively. Other parameters 
are default in QCD factorization estimations. 



small ro the quantity S{KsKs) is nonzero 

2u 



S{KsKs 



-.rr, cos{6d — 5p) sin j3. 



a/cu^ — 2a; cos 7 + 1 

For ra ~ 0.5 and — (5p ~ 0°, the CP asymmetry can reach to be S{KsKs) ~0.5. 



(24) 



III. CONSTRAINING CHARMING PENGUINS FROM GLOBAL FIT 

Let us go a step further for a model independent determination of charming penguin. 
Since the charming penguin and QCD penguin always come together, distinguishing the 
two relies on their different interference with tree type diagrams T, C and the different 
contributions for AS = and 1 modes. It also depends heavily on the precision of the 
experimental data on CP asymmetries. To isolate SU (3) breaking and possible new physics 
effects, we shall proceed in two steps: i) Fit only to AS* = modes including 11 available 
data points in tttt and KK in SU{3) limit, which is a minimal set to determine the seven 
parameters T, C, 6c, P, Sp, D and 5d- The data for K^K~ is excluded, because it constrains 



only the annihilation diagram E. 
relation from the isospin analysis 



?he subleading Pew is fixed to T + C through the SM 
The advantages to use this data set are 



that the /S.S = 0, b ^ d modes are expected to have less SU{3) breaking and less affected 
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by possible new physics. The main disadvantage is that the accuracy is limited by fewer 
data points, ii) Fits to both AS = and 1 modes which includes 19 available data points 
in TTTT, KK and nK modes. Using the approximate SU{3) symmetry, the fit accuracy is 
greatly improved. The stability of the result is checked by fit with different SU{3) breaking 
schemes. As already mentioned, the potential puzzles in i? ^ tttt and nK modes can be 
divided into hadronic dynamics related and new physics related ones. Since the current 
data show a significant reduction of the ttK puzzle in decay rates, the implication of new 
physics beyond the SM is mostly related to the low S{7!'^Ks) which remains to be confirmed 
by future experiments. We shall exclude this data point in the fits because they have little 
effects in determining the charming penguin and shall discuss it separately. 

There already exists a number of global fits to charmless B decays using fiavor diagram- 



matic methods 



36 



41 



42 



43|, 



44| and fiavor SU{3) invariant amplitudes 



431, which 



focus on using the data as an independent determination of the weak phases in the unitarity 
triangle. A recent analysis js^l shows an essential agreement with the global CKM fit results 
on the profile of UT [l9j. Since the purpose of the present work focuses on the charming 
penguin induced FSI, we simply take the values of Vub , Vcb and the weak phase from the 
global CKM fit given in Eq.(l6]) and (JT]) as inputs to further reduce the uncertainties. In fact. 



it was shown in ref. 



36| that the resulting weak phase 7 from a model independent global fit 



is consistent with the standard model and insensitive to the various cases, such as the new 
physics effects in electroweak penguin sector, the SU(3) fiavor symmetry breaking effects in 
strong phase and the charming penguin effects. As a convention, all the Brs are in units of 
10~^ and the phase angles are in gradient and arranged in the range of (— vr, +7r). 



A. Fit to TTTT and KK modes 

The fit to Tin and KK decay modes are summarized in Fit.la(b) in TabHTl For comparison 
purpose, we give in Fit. la a determination of the leading diagrams without charming penguin 
amplitude. The result is characterized by a large C/T and also a slightly large P/T compared 
with short distance QCD factorization description 

C P 

- = 0.63 ±0.09, - = 0.19 ±0.02 (Fit.la). (25) 
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modes 


Br {xlQ-^) 


acp 


S 




5.16 ±0.22 


0.38 ± 0.07 


-0.61 ±0.08 




1.31 ±0.21 


0.48l°J2 








0.06 ± 0.05 




■K+K- 


19.4 ±0.6 


-0.097 ±0.012 






9.9 ±0.6 


-0.14 ±0.11 


(0.38 ±0.19) 




23.1 ± 1.0 


0.009 ± 0.025 






12.9 ±0.6 


0.050 ±0.025 




K+K- 


0.07 ±0.12 










-0.58+o;gg 




K-K^ 


i.OD_Q 27 


'-'•-'■^-0.18 





Table I: Experimental data for charmless = and AS = 1 B decay modes 

In fit. lb the charming penguin contribution is switched on. One sees that there is a significant 
reduction of C/T ~ 0.35 from the best fitted central values while P/T is further enhanced. 
The size of D is found smaller than that of P 

C P 

— = 0.35 ± 0.16 , - = 0.36 ± 0.19, = 0.63 ± 0.62 {Fit.lb). (26) 

The best fits favor a constructive interference between C and P which is driven by the 
large decay rate of vr^vr^. The interference between P and D is largely destructive, which 
compensates the growth of P. Due to the limited degree of freedom, the inclusion of charming 
penguin leads to large uncertainties in all the fitted parameters. The Xmm curve for D given 
in Figini shows a rather fiat dependence of xLim which sets a la upper bound of ru < 1.25. 
For a meaningful determination of D, more precise data for penguin dominant nK and KK 
modes are needed. The prediction for the yet to be measured modes are 

acp{n\°) = 0.29 ± 0.48 , ^(7r°7r°) = 0.77 ± 0.58, 
acp{K^K'^) = 0.08 ± 0.52 , S{KsKs) = 0.93 ± 0.44. (27) 

The predicted acp{K'^K^) is small but S{KsKs) is very large, which follows from the best 
fitted — (5p ~ and a large = 0.62. The uncertainties in the predictions are also large. 
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parameter Fit la 



Fit lb 



Fit 2a 



Fit 2b 



T 


0.652 ± 0.036 


0.799+°?90 


0.649 ± 0.035 


0.720 ± 0.111 


C 


0.408 ± 0.052 


n 9Q9+0.115 

— U.U71 


0.467 ± 0.044 


0.416+0-082 

— U.Oby 


5n 
"G 


' '^-0.212 


-0 983+°1°? 


—1 096 ± 132 


— 1 189 lb 229 


P 


124 ± 010 


290+°-?5? 


124 ± 006 


201 ± 119 

\j a ^ \J X 1 \j a X X 


up 


-0 450+^-??I 


—0 465 ±0 115 


—0 429 ± 052 


—0 362 ±0 116 




0.013 ± 0.001 


0.013 ± 0.001 


0.013 ± 0.001 


0.013 ± 0.001 




-0.294 ± 0.103 


-0.241 ± 0.102 


-0.448 ± 0.066 


-0.416 ± 0.078 


D 


[)[jixea) 


n 1 00+0.080 


U(jixea) 


0.078 ± 0.119 


Sd 


0{fixed) 


— 407+0-172 
'J-^^' -0.512 


0{fixed) 


-0.307 ± 0.294 


xVdof 


4.4/6 


3.8/4 


15.7/12 


15.2/10 



Table II: Hadronic parameters determined from global fit to the data. Fit. la: fit to vrvr and KK 
modes in SU(3) limit without charming penguin. Fitlb: the same as Fit. la with charming penguin 
included. Fit2a: fit to vrvr, irK and KK modes in SU(3) limit without charming penguin. Fit2b: 
the same as Fit2a with charming penguin included. 

B. Fit to vrvr, KK and ttK modes 

A stronger constraint can be obtained by including the nK modes using flavor SU{3) 
symmetry. The AS = 1 modes are penguin dominant, which constrains mostly the com- 
bination P — D, and also their relative phase from direct CP asymmetries. Although the 
current vri^ data only established the direct CP asymmetry in vr+i^", nontrivial bounds for 
other modes are already obtained. The flt with vri^ in SU (3) limit is listed in TabilTl (Fit 
2). In the case of no charming penguin (Fit. 2a), one sees an even larger C/T ~ 0.7 which is 
known to be driven by the tcK CP puzzle, and the ratio P/T slightly reduced 

C P 

- = 0.72 ±0.08, - = 0.19 ±0.01 (Fit.2a). (28) 

Note that the inclusion of vri^ modes leads to a signiflcant reduction of the uncertainty. 

The fit including the charming penguin is given in Fit. 2b. Unlike the previous fits, when 
the vri^ modes are included, the ratio C/T remains large. The inclusion of charming penguin 
only leads to a slight reduction for C/T from ~0.72 to ~ 0.58, due to the lower value of r^j 
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around 0.43 with an improved precision 
C P 

- = 0.58 ± 0.14 , - = 0.28 ± 0.17, = 0.43;°:| {Fit.2b). (29) 

Compared with Fit. la, the absolute size of D is reduced from 0.18 to 0.078. With the 
improved precision, the fit result indicates stronger constraint on ro. This can be seen from 
the curve in Fig. 6. The prediction for the CP asymmetries in tt^tt^ and K^K^ modes are 

a,p(7r°7r°) = 0.53 ± 0.15 , 5(7r°7r°) = 0.73 ± 0.16, 
a,p{K^K^) = -0.04 ± 0.22 , S{KsKs) = 0.46 ± 0.44. (30) 

The reduction of S{KsKs) is also related to the reduced r^. 

To check the SU(3) breaking effects, in Tablllllwe list the fit results for two SU(3) breaking 
scheme: one is for SU(3) breaking in T diagrams only (Fit. 3), the other one is for SU(3) 
breaking for both T and C (Fit. 4). The SU(3) breaking factor is set to fx/ = 1-22 for 
= 1 modes. The obtained results show the value of r^) is quite stable 



C P 

- = 0.58 ± 0.16 , - = 0.29 ± 0.12, td = 0.34l°J^ (Fit.Sb) 
C P 

- = 0.57 ± 0.17 , - = 0.29 ± 0.13, vd = O.Slto.tl {Fit Ah) (31) 
The corresponding curves are shown in FiglHl The SU{?>) breaking scheme in Fit. 4b 



gains the lowest in a good agreement with previous analysis [37| on SU(3) breaking. The 
prediction from Fit4.b are given by 

a,p(7r°7r°) = 0.54 ± 0.25 , 5(7r°7r°) = 0.74 ± 0.22 
acp{K^K^) = -0.02 ± 0.22 , S{KsKs) = 0.38 ± 0.19 (32) 

Thus all the Fit. 1-4 favor a small acp{K'^K^) compatible with zero as a consequence of 
6£) ~ 6p but a positive S{KsKs) = 0.3 ~ 0.4. This kind of pattern is unique for the 
charming penguin contribution, which can be used to distinguish it from other contributions 
such as possible new physics from electroweak penguin sector. Note that the current data 
of S{KsKs) are not yet conclusive. The Babar and Belle collaborations report [48|, l49j 

S{KsKs) = -l^Stofstoil , C{KsKs) = -0.40 ± 0.41 ± OM{Babar) 
S{KsKs) = -0.38 ± 0.77 ± 0.08 , C{KsKs) = +0.38 ± 0.38 ± 0.05(5e//e) (33) 
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parameter 


Fit 3a 


Fit 3b 


Fit 4a 


Fit 4b 


T 

5c 


0.651 ± 0.035 

U.40Z ± U.U4o 

-1.080 ±0.135 


0.713 ± 0.092 

n 41 Q+0.079 

-1.163 ±0.199 


0.650 ± 0.035 
n A ^7 4- n n/i Q 

U.40 / It U.U40 

-1.052 ±0.127 


0.708^°:°^? 
n A^ K+0.081 

-1.121 ±0.188 


p 

r 
dp 


n 1 0/1 4- n nnv 
U.iZ4 ± u.uu/ 

— U.oOz ± U.U40 


n 1 Qi 4- n r\Q^ 
u.iyi ± u.uyo 

— U.oii it U.iiU 


n 1 9/1 4- n nnR 

U.iZ4 It U.UUD 

n Qvo 4- n n/i a 
— U.o / Z It U.U44 


U.i»b_oog7 

n QOQ 1 n 1 01 
— U.oZo =t U.iZl 


Pew 


0.013 ±0.001 


0.013 ±0.001 


0.013 ±0.001 


0.013 ±0.001 


^Pew 

n 


— U.4:Oy zt U.UDD 


n. 41 q+0.073 

u.uo / It u.uyo 


n AO'^ 4- n nfii 

— U.4Z0 It U.UDi 


u.oyo_Q 
n nfi9 4- n 1 nn 

U.UDZ It U.iUU 


5d 


0( fixed) 


-0.260 ± 0.363 


0( fixed) 


-0.295 ± 0.388 


xVdof 


15.9/12 


15.6/10 




12.7/10 



Table III: Hadronic parameters determined from global fit to tttt, txK and KK . Fit. 3a: fit without 
charming penguin. A SU(3) breaking factor fx/ f-w'^s. associated to tree diagrams. Fit3b: the same as 
Fit. 3a with charming penguin included. Fit4a: fit without charming penguin. The SU(3) breaking 
factor /i<'//7ris associated to both tree and color-suppressed tree diagrams. Fit4b: the same as 
Fit4a with charming penguin included. 

The Babar result for S{KsKs) favors a value outside physical region and has different sign 
for C{KsKs). In the Fit. 2-4, including the two free parameters D and 6d only leads to 
slight reduction of the xLin f^om 4.5 to 3.9 for Fit. lb (from 12.2 to 11.8 for Fit. 4b). Thus 
the charming penguin does not play an important role to improve the agreement with the 
data. The best fitted C/T remains large around 0.6. The charming penguin can not play a 
significant role in reducing the vrvr puzzle. 

We have checked the infiuence of the measurement of S{tt^Ks)- Including this piece of 
data leads to a big increase of the but all the best fitted parameters remain unchanged. 
For instance, in the SU(3) breaking scheme of Fit. 4b, we get = 17.7, C/T = 0.58 ± 0.12 
and = 0.34;1^q!34- As mentioned before, a low S{n^Ks) can hardly be accommodated 
within the SM, and can be a signal of new physics. A possibility is that Pew carries a 



large CP phase [35 
discussed in Refs 



53 



50 



51 



52|. Some strategies for extracting new physics parameters are 



5J|. On the other hand, the vrvr and nK CP puzzle are more relevant 
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0.05 0.1 0.15 0.2 0.25 



0.3 

|D| 



Figure 6: as functions of D. The solid, dashed, dotted and dot-dashed curves correspond to 
Fit. lb, 2b, 3b and 4b respectively. 

to the low energy hadronic dynamics and should be investigated separately. 
IV. DISCUSSIONS AND CONCLUSIONS 



An equivalent way to see the charming penguin effects is to take Ptu as a free parameter 
not equal to P and avoiding introducing the amplitude D. A global fit including Ptu was done 
a few years ago which favored a large difference between Ptu and P, hence a large charming 



penguin was implied [4]J- However, the data have been significantly changed over the years. 
The main changes in the data include i) a reduced vrvr puzzle from i^jrvr = 0.83 ± 0.23 
to the current value of 0.51 ± 0.08 ; ii) a reduced but more precise value of ^(Tr+Tr^) from 
-0.70±0.30 to -0.61±0.08; iii) a more precise a^piTf+K') from -0.09±0.03 to -0.097±0.01 
and acp(vr°_ft'~) from 0.00 ± 0.12 to 0.05 ± 0.025. The updated data are moving towards a 
much stronger constraints on the charming penguin. Our present conclusion is therefore 
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different from the previous one. 

In conclusion, we have found strong constraints to the charming penguins from its corre- 
lated contributions to B nn, nK and KK decay modes. These correlations are illustrated 
by adding the charming penguin amplitudes to these decay modes while assuming that 
other hadronic amplitudes are short-distance dominated. The charming penguin contribu- 
tion has negative correlations between Br{Tr'^7r') and Sr(7r°7r°), acp{7r'^Tr~) and acj,(7r°7r°), 
O'cpi'^'^K'') and acp(7r°X°) respectively. Positive correlations are found between acp{Tr'^K~') 
and acp{Tr~K^), acp{7i'^K^) and acp{n~K^). These correlations are unique nature of the 
charming penguin, and can be used to distinguish its contribution from the others. Using 
the latest data and assuming the approximate flavor SU (3) symmetry, the size of charming 
penguin is constrained from a global fit. The results favor a small 

td — 0.3 — 0.4 and 6d — 5p. 

which makes it unlikely as a solution to the tttt puzzle. The color-suppressed tree amplitude 
remains large C/T ~ 0.6. The time-dependent CP asymmetries in tt^tt" and KsKs modes 
are highly sensitive to the charming penguin. We have found that charming penguin leads 
to a sizable S{KsKs) ~ 0.3 while keep acp{K^K^) compatible with zero. 
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